function ans = f(x)

ans = 0.0333*x.^6 - 0.3*x.^5 - 1.3333*x.^4 + 16*x.^3 - 187.2*x + 172.9;
   
